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The active control of rotordynamic vibrations and stability by magnetic bearings and 
electromagnetic shakers has been discussed extensively in the literature. These devices, though 
effective, are usually large in volume and add significant weight to the stator. The use of 
piezoelectric pushers may provide similar degrees of effectiveness in light, compact packages. 
This paper contains analyses which extend quadratic regulator, pole placement and derivative 
feedback control methods to the “prescribed displacement” character of piezoelectric pushers. 
The structural stiffness of the pusher is also included in the theory. 

Tests are currently being conducted at NASA Lewis Research Center with piezoelectric 
pusher-based active vibration control. The paper presents results performed on the NASA test 
rig as preliminary verification of the related theory. 
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NOMENCLATURE 
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[K] 

[K d ] 

[K dd ] 

[ Kp ] 

K v 

K, 
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[M] 
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Coefficient matrix associated with {X} in state space 
Coefficient matrix associated with {i^} in state space 
Coefficient matrix associated with {U} in state space 
Damping matrix 

Uncoupled velocity feedback damping matrix 

Feedback positive active damping 

Feedback damping matrix 

Damping coef. of the piezoelectric stack 

Coef. matrix associated with the output vector 

Eccentricity 

Imbalance forces due to mass imbalance 
External forces (disturbance) 

Force produced by the i th pusher 

System gain matrix 

Identity (unity) matrix 

Polar moment of inertia 

Moment of inertia 

Minimize the performance index J 

Stiffness matrix 

Stiffness matrix including the pusher stiffness 

Pusher stiffness matrix 

Feedback stiffness matrix 

Preload spring inside the pusher 

Stiffness of the stack of piezoelectric discs 

Number of piezoelectric pushers 

Mass matrix 

Number of degrees of freedom 
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[P] 

: Matrix solution from Riccati equation 

[Q] 

: Weighting matrix associated with state vector 

r 

: Number of observer’s output 

[R] 

: Weighting matrix associated with control vector 

t 

: Number of modes used 

m 

: Control force matrix 

m 

: State vector 

{ Y } 

: Output vector 

m 

: Space coordinates 

z* 

: Pusher tip displacement 

w 

: Prescribed displacement of the pushers 

Ci 

: i th modal damping 

{0 

: Modal coordinates 

[$ d ] 

: Mode shape matrix 


: Transpose matrix of mode shape matrix 

m 

: i th natural frequency 

[diag( )] 

: Diagonal matrix 


INTRODUCTION 

An increasing amount of research is being devoted to developing effective active vibration 
control packages for rotating machinery, machine tools, large space structures, and in robotics. 
The advantages of active control over passive, i.e., absorbers and dampers, is the versatility of 
active control in adjusting to a myriad of load conditions and machinery configurations. This 
is clearly illustrated when one considers the very narrow bandwidth that a tuned spring mass 
absorber is effective in. 

Electromagnetic shakers and magnetic bearings have been used for actuators in the majority 
of the active vibration control research mentioned in the literature. Schweitzer (1985) examined 
the stability and observability of rotor bearing systems with active vibration control, and 
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presented an analysis which related force and stiffness to electrical and geometrical properties of 
electromagnetic bearings. 

Nikolajsen (1979) examined the application of magnetic dampers to a 3.2 meter simulated 
marine propulsion system. Gondholekar and Holmes (1984) suggested that electromagnetic 
bearings be employed to shift critical speeds by altering the suspension stiffness. Weise (1985) 
discussed proportional, integral, derivative (PID) control of rotor vibrations and illustrated how 
magnetic bearings could be used to balance a rotor by forcing it to spin about it’s inertial axis. 
Humphris et al. (1986) compared predicted and measured stiffness and damping coefficients for 
a magnetic journal bearing. 

Several papers describe active vibration control utilizing other types of actuators. Feng 
(1986) developed an active vibration control scheme with actuator forces resulting from varying 
bearing oil pressure. Heinzmann (1980) employed loud speaker coils linked to the shaft via ball 
bearings, to control vibrations. 

This paper develops theory and shows test results corresponding to incorporating piezoelec- 
tric pushers as actuator devices for active vibration control. The usual application for these de- 
vices is for obtaining minute position adjustments of lenses and mirrors in laser systems (Burleigh, 
1986). In the proposed application the pushers force the squirrel cage - ball bearing supports 
of a rotating shaft. The induced vibration counteracts the unbalanced vibration of the shaft by 
contributing active damping to the system. The paper presents active vibration control theory 
and test results for the piezoelectric pushers. To the authors’ knowledge this represents a new ap- 
plication of piezoelectric actuators although there has been previous applications to the bending 
vibration of non-rotating beams using 'layered piezoelectric materials (Tzou, 1987). 

THEORY 

The piezoelectric pushers consist of a stack of piezoelectric ceramic discs which are arranged 
on top of one another and connected in parallel electrically. The stack expands in response to 
an applied voltage which causes the electric field to point in the direction of polarization for 
each disc. The extension of the pusher depends on the number and thickness of the discs and 
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the force depends on the cross sectional area of the discs. Figure 1 shows a sketch of a pusher 
and the corresponding ideal model. The model consists of a prescribed displacement (a) which 
is proportional to the input voltage and a spring ( K s ) representing the stiffness of the stack of 
piezoelectric discs. The stiffness K p is a preload spring which is typically 0.001 to 0.01 times the 
stiffness K a . The figure shows that the device has a bilinear spring unless the tip is sufficient ly 
preloaded to maintain a zero gap at all times. The model utilized in the upcoming analysis 
neglects nonlinearities in the electrical and structural characteristics of the devices and damping 
(C s ) in the piezoelectric stack. 

If pushers are attached to m distinct degrees of freedom of the model, it's matrix differential 
equation may be partitioned and rearranged into the form, 


[-^f](NxiV){^}(Nx]) + [C](iVxJV){^}(Arxl) + [■^ £> ](^x^){^}(ATxl) = 

{Fd(<)}(at x i) - [K DD ](HxM){ a }(Mx 1) (1) 

where N and M are the number of degrees of freedom of the rotor and the number of piezoelectric 
pushers, respectively, and 
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The matrices [M], [C], and [K] are the mass, damping, and stiffness matrices of the rotor bearing 
system without the pushers installed, as defined in (Palazzolo, 1983). The Ki are the effective 
stiffness of the pushers, which from Figure 1 are Ki = (1 /K' a i + The stiffness Ki is 

inserted at the degree of freedom whose motion is the same as the tip motion of the corresponding 
pusher. The parameter a, is the prescribed internal displacement of pusher i, which is assumed 
to vary linearly with input voltage. 

The following portions of the paper provide the mathematical means for incorporating the 
piezoelectric pusher model into three standard active vibration control algorithms. 

Part I: Optimal Control 

Define the modal transformation 


{£}(Nxl) = [$ D ](Arxt){£}(txi) (5) 

where t is the number of modes used in the modal space and [$°] is the mode shape matrix for 
the system that includes the pusher stiffness. Substitute Eq.(5) into Eq.(l) and premultiply by 
[$-D]T, the transpose matrix of [$' D j. Furthermore, assume that the system is proportionately 
damped, and the mass matrix has been orthonormalized; 

[*°] T [Afl[§ D ] = [/](,„,) (6) 

Then, from the orthogonality conditions, Eq.(l) can be expressed as 

{{} + [<iia S (2(f n? )]{{> + [*aj([n?] J )]{«} = [$°] T {F D } - [*°] 7 '[Jf DD Ha} (7) 


Eq.(7) can be written as a first order (state space) form by adding an identity equation as 
follows. 
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Multiplying Eq.(8) by the inverse of the leading coefficient matrix yields 


{X m }( 2 |xl) — [-‘4 Dm ](2tx2t){-^ T "}(2txl) + [#F m ](2tx/V){-FD}(Nxl) 

+ [£t/ m ]( 2 txM){<*}(Mxl) 


( 9 ) 
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(ii) 


( 12 ) 


( 13 ) 


Since it is impractical to measure the displacement and velocity at all the system degrees of 
freedom (dof), an observer system may be constructed to estimate the state vector from a smaller 
number of measurements. This approach, however, may not be feasible for a large rotordynamic 
system involving many degrees of freedom and excitation at high frequencies. Another approach 
is output feedback, i.e., limit the control measurements to only those defined in the output vector, 
{T}, defined by 

{Z}\ 


{^(rxl) = [Cs](rx2N) ( ••• I 

/ (2Afxl) 

The control for the prescribed “internal” displacements of the pushers is 


(14) 


W(Mxl) = -[^(Afxrji^^rxl) 


(15) 
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where r is the number of sensors (velocities and displacements), and M the number of pushers. 
Substituting Eq.(14) in Eq.(15) and converting this equation to the modal space with Eq.(5) 
yields 

V] i [o] Wtt>' 

«> 


= -lG’][Cs] 

or in abbreviated notation 


(16) 


[0] : [* D ] 
{ ft }(Mxl) = -[G f ](Mx2t){^’}(2txl) 


where 


[G](Mx2t) = [ G'](Mxr)[Cs](rx2N ) 


l* D ] 

[ 0 ] 


[ 0 ] 

l* D l 


(17) 


(18) 


(2JVx2t) 


The objective of output feedback control in modal space can now be identified as obtaining 
the gain matrix [G’\ in Eq.(15) that suppresses the modal coordinates & in Eq.(5). The linear 
quadratic regulator problem goes one step further and also simultaneously reduces the required 
control displacements {«;}. The performance indek to be minimized is defined as 

mm{ a }J = mtn w f ({X} T [<?]{X} + {a} T [R]{a})d< ( 19 ) 

Jo 

where [Q] and [R] are symmetric, positive-definite weighting matrices which govern the relative 
importance of minimizing the modal coordinates & and the prescribed pusher displacements, 
otk- Diminishing [R] will result in larger pusher “internal” displacements but smaller vibrations 
(governed by the &). Optimal control theory produces the solution for Eq.(19) in the form of 
Eq.(17). The gain matrix [G\ is computed from 

[^](Mx2t) = ](Af x2t) [■^ > ](2tx2t) (20) 

where [P] is obtained as the solution matrix to the algebraic (steady state) Ricatti equation 
(Palazzolo, 1988). The steady state (algebraic) Ricatti equation is 


[p\[A Dm \ + {A^m - [P][Bn[jr i [BS , ”riP] + wi = [ow,,. 


iDmiTr 


)Dm 


-1 r oDmiTr 


( 21 ) 
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where the matrices [P], [ A Dm ], and [Q] are 2t by 2t, the matrix [Bfj m ] is 2t by M, and the matrix 
[R] is M by M. 

If Eq.(18) is to be solved exactly for G', the number of measured outputs r (sum of velocities 
and displacements) must equal twice the number of modes, t, used in the modal space. Define 
the matrix [0] as 

( ; 1°) ) 

[0}(rx2 1) = [Cs](rx 2 N) I I (22) 

V [ o ] : [* D ]/ ( 2 N * 2t) 

Assuming that r=2t and [/?] is nonsingular, the gain matrix [G'J can be derived from Eq.(18) and 
Eq.(20) as 

[G , ] = [R]- i [BB m } T [Pm- 1 (23) 

The prescribed “internal” displacement of the i th pusher can now be expressed from Eq.(15) and 
Eq.(23) as 

* = -£<?«»< < 24 ) 

2=1 

The force produced by the i th pusher is 

F Pi = ki(Z; - ai ) (25) 

where Z* is the pusher tip displacement. Substituting Eq.(14) and Eq.(15) into Eq.(l) yields 
[-M](iVxJV){-Z}(JVxl) + [C](JVxAT){^}(JVx 1) + [-^^(Nx^l^^Nxl) = 

/m 

{■F!dO)}(Nx 1) + [-^ I?I> ](ArxAf)[C r, ](Mxr)[C's](rx2N) I 

\{Z} 

{FD(t)}(Nxl) + [C , T’](NxN){-^}(Nxl) + [Kf](NxN){Z}(Nx1) (26) 


1 (2 Nxl) 
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[K dd ][G')[Cs] 


(27) 
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Thus, the closed loop equilibrium equation becomes 


|-M](/tfxJV){^}(Nxl) + [[C] " [£f]](JVxJV) {2}(Nx1) + 

[[Jf D ] - [if,]] {Z}(NX1) = {F D (i)}(Nx., 


( 28 ) 


This equation is very useful for conducting rotordynamic simulations with feedback control 
utilizing piezoelectric pushers. 

Part II: Pole Placement 

Similar to optimal control, Eq.(l) can be written in first order (state space) form as 
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Premultiplication by the inverse of the leading coefficient matrix yields 
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This equation is written in abbreviated notation as 

{-^}(2Nxl) = [>l £) ](2Nx2N){^}(2Nxl) + [Bp ](2NxN){Fd}(Nx1) 
+ [-B Q ](2JVxM){ Q: }(Afxl) 
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The output vector, {T}, and the output feedback control displacement, {a}, have the same 
definitions as in Eq.(14) and Eq.(15), respectively. Substituting Eq.(14) into Eq.(15) yields 

{ a }(Mxl) = ~[G , ](Mxr)[C , s](rx27V){^}(2Nxl) (33) 

Substitute Eq.(33) into Eq.(31) and rearrange 

{ x } = ([A D ] - [B„][G'][C S ]){V} + [JJFKFd} (34) 

Consider the unforced system with 

{Tb} = {0} k {X} = [*]e xt (35) 

The characteristic equation for the closed loop system becomes 

det(A[I]( 2 jVx 2 N) - [■'4 £> ](2Nx 2JV) + [#a]( 2 NxAf)[G'](Mxr)[Cs](rx 2 N)) = 0 (36) 

Assume that A is not an eigenvalue of the open loop system, then Eq.(36) can be rewritten as 

det(A[J 2 iv] — [■4 £> ])(2Nx2iV) det([/]( 2 ^ x2 N) + - [^ £> ]}( 2 JVx 2 JV) 

[■5a](2NxM)[G'](Mxr)[Cs](rx2N)) = 0 (37) 

Apply the following determinant identity (Brogan, 1974) 

det([I 2 ;v] + [PP]( 2 NxM) \QQ\(M x 27V)) = det([I M ] + [<?<3](Mx2N)[T’P](2iVxM)) (38) 

to Eq.(37) which then implies 

det([/Af] + [G ! '](Mxr)[C , s](rx2N){^[^2N] “ [^ D ]}(1 nx2N)[^“](2NxM)) = 0 (39) 

Define 

[$(^)](2Nx2A0 = {A[/ 2 n] - [^ D ]}(’ 2 Nx2N) ( 40 ) 
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and 

W^)](rxM) = [Cs](rx2N)[$(ty](2Nx2N){B a ](2NxM) (41) 

Therefore, Eq.(39) becomes 

det([/jvf] + [G'](Mxr)[^(^)](rxM)) = 0 (42) 

Eq.(42) implies that the solution will be satisfied if the column vectors of 

[Im] + [G'm A)] (43) 

are linearly dependent. Following (Fahmy, 1982) and (Stanway, 1984) this condition is expressed 

by 

{[I M ] + [G"][*(Ai)]) {/»}(Atfxi) = {0} ( Mxi) (44) 

for some {/,•} G R m . For r prescribed values of A, Eq.(44) can be expressed as 

[fl'fc* • • 'frj(Mxr) + [^[^(^l )/l ' * * - :$ W/r](rxr) = [0](AfXr) (45) 

or in abbreviated notation, 

W(Mx r) + [G'][TF]( rXr ) = [0]( rxr ) (46) 

The output feedback gain matrix is obtained from 

[G'i = -|£](Mxr)[H , l, (47) 

Note that r poles have to be assigned to compute the matrix inverse [W] -1 . Eq.(47) provides 
the gain matrix for prescribing r eigenvalues where r is the total number of sensor measurements 
including velocities and displacements. 

Physically, the gain matrix [(?'] should be a real matrix, thus for each prescribed complex 
eigenvalue A i the complex conjugate A*+i should be prescribed too. In this case, set {fi+i} = {/»} 
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or let both {/;} and {/t+i} be real vectors. This implies that ^ complex eigenvalue pairs (A Aj + i) 
may be prescribed. 

The pusher internal displacements are again obtained from Eq.(24). Furthermore, the 
feedback equivalent damping and stiffness matrices are exactly those shown in Eq.(27) and 
Eq.(28). 

Part III: Uncoupled Velocity Feedback Damper 

This is the simplest vibration control scheme in that it only involves the internal displacement 
of the pusher and the velocity of it’s tip. If the tip of the i th pusher is in constant contact with 
the k degree of freedom, the control law becomes 


on — GtfZii , i — 1 » 2, ... ,i7i 
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Substitution of Eq.(49) into Eq.(l) produces the closed loop equilibrium equation; 


[M](N X N){Z }( NX !) + + [■K’ D ](NxiV){^}(Nxl) — {^b(<)}(Arxl) (50) 


where 



\C) = [C] + 



V. 



.../ 
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and 

C* = —kiG'a, i = 1,2, ... ,m (52) 

Eq.(52) shows how positive active damping may be added into the rotor bearing system via the 
simple control law in Eq.(48). 

EXPERIMENTAL RESULTS 

An air turbine driven rotor rig was instrumented to check the Uncoupled Velocity Feedback 
Damper theory described in the previous section. The piezoelectric actuators utilized in the 
tests were Burleigh Pusher PZL-lOO’s, being driven by Burleigh PZ-150/150M Amplifier Drivers. 
Figure 2 shows a typical voltage vs. tip displacement plot for this arrangement. The curve 
in this plot provides an approximate description of the internal displacement (a) vs. voltage 
relation, since the tip is unloaded and the preload spring in Figure 1 is very light (~ 20,000 
N/M). Therefore it is assumed that the voltage sensitivity for a is 5^=-57, 000.0 V/M. This 
value was nearly constant for the three pushers that were tested. 

Load deflection characteristics of the pushers were obtained by securing each one in a solid 
cylinder, applying load to the protruding tip of the pusher, and measuring the tip deflections 
shown in Figure 3. Repeated tests with 3 separate pushers yielded an average stiffness of 
approximately 3.5 xlO 6 N/M. 

Figure 4 shows a simple sketch of the test rig, which consists of a 2.5 cm diameter shaft, 
61.0 cm in length; a 14.0 N overhung disc, 13.0 cm in diameter; and two squirrel cage mounted 
ball bearings. The outboard bearing is externally forced by an orthogonal pair of piezoelectric 
pushers, which are in turn positioned opposite to the two eddy current displacement probes d 3 
and d 4 . The control law in Eq.(48) applied to the test setup becomes 

&hor = —C*hoT ^4 (53) 

&ver = G ver d,3 (54) 
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The horizontal and vertical active damping were set equal in the control arrangement of Figure 5. 
The figure outlines how an effective damping value can be computed once the probe and actuator 
sensitivities and actuator stiffness are known. The calculations for this test setup show that the 
“active” damping coefficient is estimated to vary according to 

C A = (275) x (55) 

where G is the amplifier gain in Figure 5. The rotor was carefully balanced, and then intentionally 
unbalanced by a known amount (10.2 gm-cm), in order to compare the test results with those 
predicted by an unbalance response computer program. 

Figure 6 and Figure 7 show the test vibration amplitude vs. speed plots for the disc probes 
dj and d 2 in Figure 4. The family of curves is generated by switching amplifier gains in Figure 5 
and computing the effective damping according to Eq.(55). The computer simulation results for 
either probe d 2 or d 2 , over the same range of damping values is shown in Figure 8. A comparison 
of Figures 6, 7 and 8 show that although the test damping is less than the predicted value from 
Eq.(55) a considerable amount of damping (10,000. N sec/M, 57 lb sec/in) is still produced. 

Figure 9 and Figure 10 show the measured unbalanced response plots at the vertical (d 3 ) 
and horizontal (e^) bearing probes, respectively. Figure 11 shows the theoretical response for the 
same probes. The results again indicate that Eq.(55) overpredicts the active damping, however 
by comparing the plots the pushers do provide approximately 14,000. N sec/M (80.0 lb sec/in) 
damping at the highest amplifier-gain setting. 

The above comparison was also performed including gyroscopics in the theoretical model. 
This effect only caused minor changes in the predicted response. 

SUMMARY AND CONCLUSIONS 

This paper has examined the possible use of piezoelectric pushers for active control of 
rotor-bearing system vibrations. Although their most common application is currently micro- 
positioning of laser system mirrors and lenses; their stroke, force and frequency response make 
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them potentially very useful for vibration control. The obvious advantage of piezoelectric pushers 
over other actuators is their compact size and light weight. 

Three active theories were extended to treat the special prescribed “internal” displacement 
(a) character of piezoelectric pushers. The Uncoupled Velocity-Feedback Damper theory was 
then tested by comparison to experimental results. The study showed that although the theory 
overpredicts the amount of damping produced by the pusher the actual level is significant, being in 
the range of 10,000-14,000. N sec/M (57.0-80.0 lb sec/in). The discrepancy between predicted and 
measured results most likely arises from the nonlinearity and hysteresis in the voltage-deflection 
and load-deflection characteristics of the pushers, and from the neglect of the structural damping 
of the piezoelectric stroke, by the theory. We are currently working with the pusher manufacturer, 
Burleigh Inc., to produce pushers that have reduced nonlinearity, hysteresis and stack damping. 
Construction of a feedback box to perform optimal control (OC) and pole placement control 
(PPC) is in progress. The device will accomodate up to 12 sensors and produce outputs to two 
actuators, which corresponds to a 2x12 dimension for the [(?'] matrix in Eq.(15). Test results 
from application of OC and PPC will be forthcoming in the literature. 
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Figure 1 Sketch of piezoelectric pusher and corresponding analytical model. 
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INPUT VOLTAGE (V) 



Figure 2 Typical pusher tip deflection vs. input voltage. 
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Figure 3 Typical pusher tip deflection vs. pusher tip load. 
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Figure 4 Diagram test rig with piezoelectric pushers. 
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Sd — -3t? / mil = ll,8l0.v/m 


Prox. Prox. 




8.6 x 10“ 4 d/dt 


J Kd \ wrticd pie*o. 


Amplifier 

Gain=G 


Amplifier 

Gain =G 




Synchronous 
Tracking Filter 


(Out =0.44 x (input) 


Synchronous 
Tracking Filter 


■ 

I 


Pusher 


Driver 


Uncoupled Velocity Feedback; a = G* d. 


Active Damping 


C A = K a <? 


Actuator input voltage for a displacement, V Al is 

V A = S A aj 

FVom (a) and substituting the data shown in Figure; 


V A = 5 A CTd = (8.6 x 10“ 4 d)G(0.44) 


Hence 


S A £±d = S D x G(3.8 x 10- 4 )d 

AX 


= 3.8 x 10” 4 x C x Jif^ x *ec 

Sa 


Ca = (3.8 x 10 ~ 4 sec) ~~~~r~(3.5 x 10 *N/m) x G 
/ 57,000t>/m v 

___ -Nsec Absec. 

= 275 x G (1.57 x < 7 — — ) 

m tn 


Figure 5 Outlines of system setup and the related equipment coefficients. 


107 











RESPONSE (MM) (P-P) RESPONSE (MM) (P-P) 


PRESCRIBED PUSHER DISPLACEMENT 
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Figure 6 Prescribed pusher displacement at probe d\ (DOF 1). 
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PRESCRIBED PUSHER DISPLACEMENT 

D1 PROBE (DOF(1 )) 



Figure 8 Computer simulation results of prescribed pusher displacement, d t 
probe. 
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PRESCRIBED PUSHER DISPLACEMENT 

D4 PROBE (DOF(7)) EXPERIMENT DATA 



Figure 10 Measured unbalanced response plot at horizontal, d,\, probe. 
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